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Abstract 

We derive model-independent constraints on four-fermion contact interaction-type dynam- 
ics from the published preliminary LEP2 experimental data on e + e~ annihilation into 
fi + fi~ and t + t~ pairs, measured at different energies between 130 and 207 GeV. The basic 
observables are chosen to be the total cross section and the forward-backward asymme- 
try, and the analysis realistically takes into account data uncertainties and correlations 
among measurements at the various energies. The combination of data from different en- 
ergy points plays an important role in the determination of regions allowed for the contact 
interaction coupling constants. In contrast to the more common one-parameter analyses, 
we only obtain constraints on pairs of parameters rather than limits on individual ones. 
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1 Introduction 



Many Standard Model extensions envisage a dynamics acting at one (or more) large mass 
scales A 3> Mw, such that the relevant states exchanged among quarks and leptons, 
having a mass proportional to A, are so heavy that they could not be directly produced at 
accelerator energies. The most familiar case is represented by quark and lepton composite 
models PjJ [2j, but there are numerous other examples. However, such new interactions 
could manifest themselves by indirect, virtual, effects represented by deviations of the 
measured observables from the Standard Model (SM) predictions. If some deviations were 
effectively observed experimentally to a given significance level, one could try to derive from 
the data numerical information on the parameters (masses and coupling constants) of the 
non-standard models and, eventually, to select the viable one. In the case where, instead, 
no deviation from the SM predictions were observed within the experimental accuracy, one 
can set numerical bounds and/or constraints on the parameters characterizing the new 
interactions and, in particular, on the relevant mass scales A. This information should also 
be of phenomenological interest, in the exploration of non-standard interactions. 

In the spirit of "effective" theories, exchanges of very heavy objects in reactions of 
quarks and leptons can be parameterized by a contact interaction, representing the "low 
energy" expansion of the transition amplitude to leading order in the small ratio \fs/A 
(\fs being the cm. energy). The explicit form of such contact interaction Lagrangian (CI) 
depends on the particles participating in the reaction under consideration. Specifically, we 
consider here the electron-positron annihilation: 



with f = fj, and r, and the relevant precision data at LEP2 for 130 < y/s < 207 GeV, 
published in Ref . j2] , where the results of the four experimental collaborations are combined. 
Such high precision data can be regarded as a powerful tool to severely test manifestations 
of non-standard interactions through deviations from the SM predictions. In particular, 
we are interested in deriving, from those data, constraints on the e e f f contact-interaction 
Lagrangian [2]: 



where a, (5 = L, R denote left- or right-handed fermion helicities, and the parameters 
e Q/ 3 specify the chiral structure of individual interactions and determine the size of the 
deviations from the SM predictions. One can introduce the previously mentioned large 
mass scales by \e a p\ = l/A^g, and coventionally fixing g^/^n = 1 as a reminder that, 
as a compositeness remnant force, this interaction would become strong at y/s ~ A a p. 
However, as remarked above, more generally the scales A a p define a standard to compare 
the sensitivity of measurements to the various kinds of new interactions, see, e.g., [UGI]. 

In practice, the situation is complicated by the fact that, for a given fermion flavor 
/, Eq. (J2J defines four individual and independent models (basically, the combinations of 
the four chiralities a, (3 through the e's) and, in principle, the general contact interaction 
could be any linear combination of these models. Thus, the aforementioned deviations of 
the cross section from the SM predictions may simultaneously depend on all four-fermion 



e + + e 



/ + /, 





(2) 



2 



effective couplings and, if only one value of the cm. energy were available, the straightfor- 
ward comparison of deviations and experimental uncertainty could produce only numerical 
correlations among the the different CI couplings, rather than separate, and restricted, al- 
lowed regions for these parameters in the parameter space around the SM limit e a p = 0. 
Moreover, negative interference of CI and SM amplitudes in the cross section might con- 
siderably weaken the bounds. 

The simplest and commonly adopted procedure consists in assuming non-zero values 
for just one of the e a p at a time, and in constraining it to a finite interval by essentially a 
X 2 fit analysis of the measured cross sections and forward-backward asymmetries, while all 
the other parameters are set equal to zero [IHOj- In this way, only tests of specific models 
can be performed. 

On the other hand, it would be desirable to perform a more general kind of analysis of 
the data, that simultaneously includes all terms of Eq. (J2J) as free, potentially non- vanishing 
independent parameters and, at the same time, allows to disentangle their contributions 
to the basic observables in order to derive separate constraints within finite regions around 
the SM limit. 

In cases where only one value for the cm. energy is available, such as for the planned 
e + e~ Linear Collider JTj. a solution is represented by the initial electron beam's longitudinal 
polarization, that would enable to experimentally extract the individual helicity amplitudes 
of process (JTJ), by definition directly related to the individual e e / /contact couplings e a p 

pin!. 

Such a procedure cannot be applied to the data from LEP, with unpolarized electron 
and positron beams. However, in this case, the cross sections of processes (JQ are measured 
at LEP2 over a range of a/s values wide enough that the energy dependence of the devia- 
tions, entirely determined by well-known SM parameters, can be exploited to restrict the 
bounds to limited regions in the CI parameter space, and in this way to perform an analysis 
of the new interaction, model-independent in the sense indicated above. This observation 
was used for a global analysis of data at the energies of LEP1, LEP2 and TRISTAN in 
Ref. (lOl- The analysis presented here uses exclusively the most recent higher statistics 
LEP2 data, combines the two channels /i + fi~ and r + r~ and the results of the four exper- 
iments, and accounts for, among other things, the correlations among the measurements 
at the different energy points. The basic observables will be the "conventional" ones, 
namely, the integrated cross section er(s) and the forward-backward asymmetry Apb(s), 
whose experimental values are tabulated in Ref. [3]. 

Specifically, in Sec. 2 we will give the basic definition of helicity amplitudes and the 
formulae relevant to o and Afb for the processes of interest here, and in Sec. 3 we shall 
present the model-independent analysis of LEP2 data and the resulting constraints on CI 
couplings. Finally, Sec. 4 will be devoted to some concluding remarks and an application 
of the method to a model example. 

2 Cross section and helicity amplitudes 

Limiting ourselves to the cases / = /i, r and neglecting all fermion masses with respect 
to \/s, and taking into account the Born 7 and Z exchanges in the s channel plus the 
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contact-interaction term (J2J, the differential cross section of process (0) reads [TTj: 



A Q 

° [(1 + COS 0) 2 <7+ + (1 - COS 0) 2 <7_] , (3) 



dcosfl 8 

where is the angle between the incoming electron and the outgoing fermion in the cm. 
frame. In terms of helicity cross sections, o a p with a, (5 = L, R: 

°~+ = T (^LL + 0R R ) , (4) 

o-- = 7 (^lr + ctrl) • (5) 



In Eqs. (JD and (JSJ): 



O-a/3 = apt|jM Q( g| , (6) 

where er pt = cr(e + e~ — > 7* — ► Z + /~) = 47ia^ m /3s (for quark-antiquark production a color 



factor iVc — 3(1 + a s /7r) would be needed). The helicity amplitudes M. a p can be written 

as 

f s 

M a /3 = QeQf + 9a9pXz H e a /3, (7) 

C^e.m. 

where: = s /( s ~ + ^M^T^) is the Z propagator; g[ = (l( L — Q fS^r) / Sw c w and 
9r = ~Qf s w/ s wCw are the SM left- and right-handed fermion couplings of the Z with 
Syy — 1 — Cyy = sin 2 9w\ Qe = Qf = — 1 are the fermion electric charges. 
The measured observables a and Afb are given by the relations: 

d cos 9 = - [((Jll + orr) + («"lr + Ctrl)] ; (8) 

_i a cos v 4 



and 



f®\ dcr 3 
J ~ J ) dcos6> d C ° S ^ = 16 ^°" LL + CTRR ^ _ ^ LR + aRL ^ ' ^ 



Finally, their relation to cr± is given by 

a ± = °- (l±^ FB ). (10) 

Taking Eq. (JJJ) into account, Eqs. (jSJ) and © show that a and <tfb (or ^4fb) simultane- 
ously depend on all four contact interaction couplings, and therefore by themselves do not 
allow a model-independent analysis, but only the simplified one-parameter fit of individual 
models. However, a and <tfb depend on the two combinations of helicity cross sections, 
(o"ll + 0rr) and (ctlr + 0rl). Accordingly, a combined analysis of a and <tfb enables to 
separately constrain the pairs of parameters (cll;^rr) and (clr^rl)- Moreover, the com- 
bination of experimental data on a and o"fb at different values of the cm. energy allows 
to further restrict such separate bounds in a model-independent way. 

To clarify this statement and intuitively show by a simplified example the role of the 
different energy points in improving the constraints, assuming that no deviation from the 
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SM is observed within the experimental accuracies, constraints on the contact interaction 
couplings € a p can be derived from the system of two inequalities: 

\a SM+cl - a SM \ < 5a, (11) 



|^ +CI -^IbI<^fb, (12) 

where 8a and 5Ap R represent the experimental uncertainties on these observables. Taking 
Eqs. (JIJ) and (jHJ) into account, the deviations from the SM predictions in the left-hand sides 
of Eqs. (fTTjl and (fT2|) can be written as: 

^SM+CI _ ^SM = 1 [(A(Tll + A(Jrr) + (A(Jlr + A(Jrl)] ^ (13) 



^.SM+CI a 



SM 



FB FB - 16(J SM 



1 - ^A^\ (Aa LL + Aa RR ) - (l+ ~A™ ) (A<t lr + Aa RL ) 



(14) 

where Aa a , = a™ +cl - o™. 

From Eqs. pijl — ()14j) one can obtain constraints on the e a p. Specifically, the areas 
allowed to the values of the parameters are enclosed by concentric circles in the planes 
(£ll,£rr) and (€lr,£rl)- For example, the domain allowed to the pair (cll> £rr) is delim- 
ited by the circular contours: 



(** + ^M™) 2 + (e RR + ^AC) 2 = Rl (15) 
where 

Rl = (?^M™) 2 + (^AC) 2 ± -\ (16) 



and 



K 2 =(^) 2 —8a+. (17) 



a 



P t 



In the right-hand side of Eq. ([17)1 . 5a + must be expressed in terms of the experimental un- 
certainties 5a and SA-p-Q, 1 and a color factor 1/Nc is needed in the case of quark-antiquark 
production. These relations show that both the centre and the radii of the circles R± are 
determined by the values of the SM helicity amplitudes and depend on energy, while the 
width of the allowed area is determined by the experimental uncertainty of the observables. 
Therefore, in principle, the combination of two (or more) such allowed regions correspond- 
ing to different energies can lead to a reduction of the allowed region and, ultimately, to 
model-independent bounds on the contact interaction coupling constants. 

It should be stressed that, while the observables are given by two sums of helicity cross 
sections, ctll + 0rr and <tlr + orx, it does not follow that one can only obtain constraints 
on the corresponding sums of parameters, £ll 4" £rr and 6lr + crl- There is indeed a small, 
but finite, sensitivity to the individual parameters. This is due to effects proportional to 
squares of the e parameters, together with the small difference between left- and right- 
handed couplings of the Standard Model, \gl\ ^ \g R \. 

1 While <r + and <r^, as shown in Eqs. and J5J), are the most natural observables, we use instead a 
and A-fb for our analysis, since the data are tabulated (with errors) for these observables. 
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3 Data fitting and derivation of constraints 



Recently, the // Subgroup of LEPEWWG presented preliminary combined results of mea- 
surements of the four LEP collaborations using experimental data from the full LEP2 
available data set at energies from 130 GeV up to 207 GeV for the annihilation processes 
e + e~ — > ff [3]. In particular, for lepton final states / = \i and r, the set of the average 
cross sections a^, a TT and forward-backward asymmetries A^, and their experimen- 
tal errors have been given for the 12 energy points listed in Tabled 



Table 1 : Approximate average integrated luminosity per experiment and nominal centre-of 
mass energies collected during LEP2 operations [Sj. 



Ecu (GeV) 


130 


136 


161 


172 


183 


189 


192 


196 


200 


202 


205 


207 


Ant [pb" 1 ] 


3 


3 


10 


10 


50 


170 


30 


80 


80 


40 


80 


140 



The data fitting procedure used is based on the method of least squares. We introduce 
a x 2 function, which may be written in the following matrix form 

x \e) = (C> LEP2 - TH (e)) T V- 1 (e? LEP2 - C? TH (e)), (18) 

where e = (6rr, cll, crl, £lr) is the vector of C.I. parameters; O hEP2 is the vector of values 
of observables measured at LEP2 and th is the vector of their theoretical predictions; 
finally, V is the covariance matrix of the experimental uncertainties. 

The chosen set of observables, represented by the vector £} LEP2 ; contains 48 elements 
(two kinds of observable for two flavour channels and twelve energy points). The cor- 
responding theoretical predictions, th , which depend on the CI parameters e and on 
radiative corrections via improved Born SM amplitudes ^21 EH] , have been evaluated with 
^top = 175 GeV and m# = 150 GeV. Initial- and final-state radiation are taken into ac- 
count by the program ZFITTER adapted to the present case of contact interactions. 
The radiative corrections were applied using definition "2" in Ref. [3]: namely, for dilepton 
events y/s* is taken to be the bare invariant mass of the outgoing dilepton pair (as opposed 
to that of the s-channel propagator), the ISR-FSR photon interference is included and the 
signal is defined by the kinematical cut \fs' > 0.85i/s. We note that the improved Born 
amplitudes leave the form of the previous equations for the cross sections, Eqs. (J3j) and 
(fIU)l unaltered. 

As regards the 48 x 48 symmetric covariance matrix V, the diagonal entries are the 
experimental uncertainties on the observables, while the off-diagonal entries define the 
correlations between the observables as well as among the different energy points 

The least-square confidence region is determined by the condition 

X 2 (e)<xLn + XcL, (19) 

where Xmm * s the minimum value of the function x 2 ( e ) an d Xcl = 9-49 for 95 % CL and four 
degrees of freedom. The procedure of minimization x 2 ( e ) is performed using the program 
package MINUIT pS]. 
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Table 2: Central value e°, global limits (allowed intervals) obtained as projections of the 
95% CL four-dimensional region on the axes and 95% CL one-dimensional mo del- dependent 
constraints on the CI parameters. 



Parameter 
[TeV~ 2 ] 


Model independent 


Model dependent 


central value 


global limits 




0.0085 


(-0.175, 0.095) 


-0 47+ - 0071 

u - uu ^' -0.0071 




-0.0195 


(-0.187, 0.111) 


-0 00S2+ - 0078 

U.UU<JZ_ .Q 07 g 




0.0120 


(-0.225, 0.060) 


-0 0012+ 00111 
u.uuiz_ 0116 




-0.0160 


(-0.225, 0.060) 


-0 0012+ ' 0111 
u.uuiz_ 0116 



Combining the [i and r data, we show in Table |2] the components of the central value 
e° (over-all minimum of x 2 ) and the global limits (intervals (e m ; n , e max )) obtained as pro- 
jections of the confidence region on the corresponding axes. These intervals should be 
considered as global, model- independent, constraints on the CI parameters e a p. The x 2 i* 1 
the model-independent fits amounted to 41.3, for = 48 — 4 = 44 degrees of freedom: 
the probability of this result is p = 0.411 [TB]. For comparison, we give the 95% CL one- 
parameter constraints on e a p parameters for the LL, RR, LR and RL contact-interaction 
models. 

In Figs. ^121 we show the contours which bound the regions found as "projections" 
of the four-dimensional confidence hypervolume determined by (|19p on four of the two- 
dimensional planes (LL-RR), (LR-RL), (LL-LR), (LL-RL). The contours have been pro- 
duced as the line connecting all points of the plane where \ 2 takes the value Xmin + Xcl 
after minimization on the two remaining free parameters. 

These figures show obvious symmetries. First of all, in Fig. Q where we display the 
allowed regions in the (eix, £rr) and (eLR, 6rl) planes, there is an approximate "reflection 
symmetry" between 6ll £rr as well as between £lr <-> €rl. As discussed in Sec. 2, 
the observables depend on <tll + 0"rr and ctlr + ctrl, and to lowest order in the e's, this 
translates into a dependence on €ll + £rr and €lr + crl- Thus, in this approximation, 
the allowed regions would be bands at fixed €ll + crr and €lr + crl, representing strong 
correlations between pairs of parameters. The contributions of the second-order effects (in 
the e a p) delimit and curve these bands. In the case of the €lr-6rl exclusion region, the 
radius of curvature, given by an expression analogous to Eq. I|16|L is smaller than that 
of the 6ll~£rr exclusion region. This stronger bending originates from the destructive 
vs. constructive interference between photon- and Z-exchange: above the Z resonance, 
gbkxz < 0, whereas {g\f X z « (gi?Xz > 0, so \Mfg\ = \M S ™\ < \M%\ w \M$&\. 

In Fig. El we show the analogous allowed regions in the (eix, ^lr) and (£ll> £rl) planes. 
The allowed regions in the (€rr, £lr) and (crr, £rl) planes are very similar to those of this 
Fig. |2l and hence not shown. 

In the figures, the constraints on the one-parameter models LL, RR LR and RL (see 
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Figure 1: Allowed areas at 95% CL. on leptonic contact interaction parameters in the 
planes (ell, £rr) and (6lr, c R l), obtained as projections of the four-dimensional confidence 
hypervolume on the relevant plane after minimization in the two remaining parameters. 
The bars correspond to one-dimensional model-dependent constraints as discussed in the 
text. The circles correspond to the central values (see Table |2J). 



Table |2]) are represented by bars. These correspond to one- dimensional mo del- dependent 
constraints at 95 % CL. with Xcl = 3.84, obtained by varying only one parameter at a 
time with the remaining three set equal zero. The results in this case are in full agreement 
with those obtained in Ref. 

Figure El is rather different from Fig. Q but the two panels are very similar among 
themselves. This is due to the symmetric inputs, <tll + o"rr and o"lr + &rl, together with 
the fact that the linear approximation, determined by the interference between SM and 
CI couplings, provides a first, rough description of the bounds. Also, we note that there is 
little correlation between these pairs of parameters, and the allowed regions are simply 

[e aP < 5'} n [e^p < 5'} n [e 2 a/3 + e 2 a , p , < 5% (20) 

where 

5' ~ max(e LL ) or max(e LR ) or max(e RL ), (21) 
as determined from Fig. ^ Similarly, 

5 2 = [max(e LL )] 2 + {[max(e LR )] 2 or [max(e RL )] 2 }, (22) 

respectively, for the two panels. This simple shape is thus due to the lack of correlations 
among the parameters shown in Fig. 121 
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-0.2 -0.1 0.1 0.2 „ -0.2 -0.1 0.1 0.2 „ 

£ LL (TeV" 2 ) £ LL (TeV" 2 ) 

Figure 2: The same as Fig. 1 for the (cllj ^lr) and (cllj ^rl) planes. 



4 Discussion 

Our most important result is that if one does not restrict the analysis to individual models, 
the bounds on the e's are rather loose. In fact, any set of three of them (but not all four, as 
is seen from the correlations in Fig. EJ can be of the order of 0.2 TeV -2 . This corresponds 
to a scale A ~ 2.2 TeV. 

In the case of €ll and £rr, as discussed above, the orientation of the "banana" in Fig. [I] 
implies that e LL and e RR should roughly add to zero. Rather large deviations from the 
SM are allowed, provided these parameters have opposite signs. It should be noted that, 
if one assumes universality, such opposite signs can never arise from the low-energy limit 
of a vector-particle exchange, irrespective of the chiralities of the couplings. This may in 
part explain why the present bounds are much looser than those of the mo del- dependent 
analyses. 

Also, it should be stressed that we do not assume full lepton universality in this analysis. 
The muon and tau data are combined, but the couplings to those currents are not taken to 
be the same as the couplings to the electron of the initial state. If lepton universality were 
imposed, one would have the additional constraint €lr = 6rl- Also, full lepton universality 
would imply the product €ll £rr > 0, and much of the allowed part of Fig. [T] would be 
excluded. We note that there are models without flavour universality, where €ll and £rr 
can have opposite signs (see, e.g. [TTf ITS]). 

Up to this point, we have combined the muon and tau data. It is also interesting to 
study these two data sets separately. In Fig. El we show the contours which bound the 
allowed regions in the (eLL, £rr) and (6lr, 6rl) planes for e + e~ — > t + t~ (i.e., without using 
the muon data). The general shapes of these allowed regions are rather similar to those 
obtained from the combined data, but they are significantly larger. The corresponding 
allowed intervals of the e parameters are given in Table El the analogue of Table [21 The 
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model-independent global limits, for example, are looser than the combined muon and tau 
analysis by up to 40%. The muon data alone give shapes and allowed intervals quite similar 
to those in Fig. El and Table HU respectively, but narrower, essentially reflecting the larger 
total error in the r case. 



Table 3: Similar to Table El for final-state r pairs only. 



Parameter 
[TeV^ 2 ] 


Model independent 


Model dependent 


central value 


global limits 


£ll 


0.0005 


(-0.249, 0.113) 


-0 0032+ - 0118 


crr 


-0.0125 


(-0.258, 0.136) 


-0 0035+ ' 0129 
u.uuoo_ 0131 


^lr 


-0.016 


(-0.273, 0.066) 


—0 00^fi +0 - 0188 
u.uuoo_ 0202 




0.0085 


(-0.273, 0.066) 


-0 0036+ ' 0188 



In specific models, there are often constraints on these deviations e Q( g. For example, 
the Z' couplings of E e models lead to the constraints [T9^ : 

Z' x : e LR = e RL <0, e LL = 9e RR <0, (23) 

Z'^ : e LR = e RL > 0, e LL = e RR < 0, (24) 

Z' v : e LR = e RL < 0, e LL = — e RR < 0. (25) 
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These signs are given by the signs of the couplings, together with the low-energy limit of 
the propagator. The leptonic data studied here, lead to the M z > bounds: Z'\ 600 GeV: 

A, 

ZL: 330 GeV; Z' r] : 340 GeV. The corresponding bounds from all data [3J are 670 GeV, 
480 GeV and 430 GeV, respectively. 

Also, in the case of models with TeV-scale extra dimensions (with Kaluza-Klein excita- 
tions of the photon and the Z), there are relations among the couplings [20]: €lr = 6rl < 0, 
and 6ll = ^kr/^ s w — e RR < 0- For one ex tra dimension, the bound on the compactifica- 
tion scale £1] is M c > 2.2 TeV. 



>0.04 



-0.02 



-0.04 









LEP P 


'eli mint 


3ry 


























































V , 





-0.04 -0.02 0.02 0.04„ 

£rl (TeV ) 



Figure 4: Similar to Fig. ^ (right panel), but magnification of the region of small £lr and 
e^L, for Xcl = 5.99. The diagonal line corresponds to the constraint of the anomalous 
gauge coupling model. 



As anticipated in the Introduction, we shall here consider an example application, 
namely the effects of anomalous gauge couplings [T7j in the process (JTJ). We note that this 
model, which assumes universality, is characterized by two parameters, foB and Jdw- The 
deviations (jJJ) will take the form 



Idw . 2fDB \ 8Jdb 4/ 



'j DB / n „ N 

CLL — tte.m. | -pr~2 1 2 I ' 6rr ~ a e.m. 2 ' 6lr — 6rl — a e.m. 2 ' V /D / 

ZS W C W J C W C W 

where foB and fow are related to foB and fow of ref. [Ej by / = / /mf. In this model, one 
has €lr = crx, so any deviation would be restricted to lie along the dashed line in Fig. HI 
which shows a magnification of the allowed band in Fig. ^ (right panel) for Xcl = 5.99, 
corresponding to two parameters. The intersections with the allowed bounds allow us to set 
a limit on \ fr>s\ < 0.21 TeV -2 . This also amounts to a bound on €rr. From the analogue of 
Fig. H] (left panel), corresponding to two parameters, one can extract bounds on eix- Using 
Eq. (J2HJ), these can then be converted to the bounds: —1.7 TeV -2 < fow < 1-1 TeV -2 . 
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We have not analyzed the quark data, which is of poorer quality, due to the limited 
efficiency of 6-tagging, together with the problem of distinguishing b from b jets (see, 
however, Ref. [22]). 

At the Linear Collider, where polarization would be available, more observables can be 
studied, such as Alr and Alr,fb- Thus, and because of the higher energy, dramatically 
tighter constraints are foreseen [§]. 
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